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^= (54) Title: A METHOD FOR ARBITRARY TWO-DIMENSIONAL SCALING OF PHONOCARDIOGRAPHIC SIGNALS 




(57) Abstract: A heart signal comprising both first and second 
heart sounds and murmurs occurring due to various heart condi- 
tions is modeled by a sum of sinusoids automatically selected to 
represent the heart sound without a noise component. In an anal- 
ysis step the signal parameters are measured, in an interpolation 
and transformation part fi(t) and A(t) are interpolated from one 
window position to the next window position, and a syntheses 
part reconstracts the transformed heart signal. The transforma- 
tions are transformations in either time scale (with unchanged 
frequency content ) or frequency scale ( with unchanged time 
progression ), or both. Furthermore a dynamic transformation of 
an arbitrary heart signal to a fixed number of beats per minute is 
made possible. 
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A method for arbitrary two-dimensional scaling of phonocardiographic signals. 

5 The invention relates to a method for assisting medically trained persons in 

diagnosing and classifying medical conditions of the heart by means of auscultation. 
It also relates to a method and apparatus for transforming the information content of 
a phonocardiographic signal to a frequency range and duration commensurate with a 
communication channel having properties that fall outside the normal frequency 

10 range of phonocardiographic signals. 

Signals obtained by means of a transducer are phonocardiographic representations of 
sounds traditionally listened to by means of a stethoscope. Training in auscultation 
takes a long time and requires an aptitude for recognising and classifying aural cues, 
15 frequently in a noisy environment. 20-30 different conditions may need to be 

differentiated, and within each, the severity evaluated. Furthermore, there may be 
combinations among these. These factors contribute to explaining why not all 
physicians perform equally well when diagnosing heart conditions, and why it may 
be time-consuming. 

20 

Use of modem signal processing has been proposed to enhance the relevant aural 
cues, and in some cases it has been demonstrated that a transposition into a different 
frequency range or repetition of a repetitive signal at a different repetition frequency 
from that originally recorded may be advantageous. Among the patents that describe 
25 this kind of signal processing the following may be mentioned: 

US 4,220, 160 uses the heart sound picked up by a microphone to modulate a higher 
and more audible carrier frequency. US 4,424,815 is a very ftmdamental patent 
employing a range of solutions, from magnetic drum (i.e. a non-digital procedure) to 
30 a tapped delay line for stretching time while conserving frequency content. US 

4,528,689 describes a time stretching digital stethoscope that records several cycles, 
and obtains a "slow-down factor" proportional to the number of cycles recorded. 
Each part of a sound between zero crossings is repeated with opposite polarity, 
thereby obtaining a signal that sounds like original signals but lasting longer. Only 
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whole-number &ctors are possible. US 4,792,145 desoibes how frequency scaling is 
obtained by multiplying each Fourier spectral component by a factor. The original 
relationship between frequency components is retained, and it is purported that the 
audibility is improved. Similarly, the time may be scaled. 

All of the above solutions have severe limitations, either because of an increase in 
non-linear distortion and noise, or because the factor for stretching can only be 
whole-numbered. 



10 According to the invention, these disadvantages may be avoided by synthesizing a 
"clean" signal consisting of a sum of sinusoids and to perform a transformation on 
these sinusoids in order that the relative amplitudes are maintained while the 
frequency range or time axis are changed in a continuous range of scales. In a frirther 
refinement of this method, the scale may be linked to other phenomena by means of 

15 an auto-scale frmction. 



An embodiment is particular in that that the signals are converted on a running basis 
to a sinusoidal model representation, that a time and/or frequency axis scaling is 
defined and used to control the amplitudes and phases of the sinusoids, which are 
20 subsequently added to create a time and/or frequency scaled representation of the 
phonocardiographic signal. The basis of the invention is the method in which a 
correct amplitude and phase adjustment of the contributing sine generators is 
obtained. 



25 An advantageous embodiment of the invention is particular in that it comprises the 
steps of obtaining the frequency content of the signal by applying a Short Time 
Fourier Transform with overlapping segments, a Discrete Fourier Transform being 
performed on each segment, performing a frequency peak search on each segment by 
consecutive removal of the spectral components having the highest energy content, 

30 identifying each peak by its frequency value repeating the peak search, until a 

maximum number of peaks have been identified or until the energy content of the 
last peak is below a preset minimum, establishing a segment-by-segment map of 
spectral peaks, said peaks forming a track over time, optionally subjecting the 
frequency values of the spectral peaks to a multiplication, defining a synthesis frame 
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of time, based on said segments, optionally subjecting each frame to a multiplication 
of the time scale, adjusting the phase of sine generators centered on the frequencies 
of the tracks, adjusting the amplitudes of said sine generators, and summing the 
outputs of all sine generators active at any one instant for a given frame length T. 

5 

A more general embodiment is particular in that it includes the following steps: 
obtaining the frequency content of the signal by applying a Short Time Fourier 
Transform with overlapping segments, a Discrete Fourier Transform being 
performed on each segment, performing a frequency peak search on each segment by 

10 repeated removal of the highest spectral "hills" identifying each peak by its 

frequency value, repeating the peak search until a maximum number of peaks have 
been identified or until the peak level of the last peak is below a preset minimum, 
establishing a segment-by-segment map of spectral peaks, said peaks forming a track 
over time, adjusting the phase of sine generators centered on the frequencies 

15 represented by the tracks, summing the outputs of all sine generators active at any 
one instant for a given frame length T, creating a continuous output signal by joining 
consecutive frames. In the present embodiment use is made of the fact that the 
estimate of the importance of a peak (a "hill") may be obtained by other means than 
an energy measure, i.e. a different sorting criterion. 

20 

A fiirther advantageous embodiment for scaling a phonocardiographic signal on the 
frequency axis is particular in that the spectral peaks are multiplied by a fector q, 

A ftirther advantageous embodiment for scaling a phonocardio^phic signal on the 
25 time axis by a desired &ctor is particular in that the frame length is multiplied by a 
factory. 

A further advantageous embodiment for autoscaling a phonocardiographic signal on 
the time axis is particular in that the scaling factor p is set such that the frame length 
30 multiplying factor is equal to the heart rate divided by 60. This means that a stability 
of the signal will be obtained irrespective of changes in the heart rate. 
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An advantageous embodiment of the invention that ensures sufficient resolution and 
precision in synthesis is particular in that the number of sine generators is maximum 
50 in any one frame. 

5 The invention also relates to an apparatus for performing the method, being 

particular in comprising means for windowing the time fiinction, short time Fourier 
spectrum analysis means, means for searching and classifying spectral peaks, means 
for comparing phases of signals corresponding to said spectral peaks, means for 
controlling the phases of sine generators providing signals corresponding to said 
10 spectral peaks, means for controlling the amplitudes of said sine generators, and 

means for summing the signals of said controlled sine generators in order to obtain a 
synthesized and essentially noise free output signal representative of said time 
fiinction. 

15 The invention will be more fully described in the following with reference to the 
drawing, in which 

Fig. 1 shows a block diagram representative of the operations performed on an input 
signal in order that a sinusoidal model is obtained for the signal, which is 
20 summed to obtain an output signal. 

Fig. 2 shows the procedure of block segmentation. 

Fig. 3 shows the estimation of a finite number of spectral peaks. 

Fig. 4 shows results of the peaks search 

Fig. 5 shows the matching procedure 
25 Fig. 6 shows the final results of the peak matching 

Fig..? shows boundary, condition for the interpolationJiinctions, 

Fig. 8 shows interpolation of instantaneous phase function and amplitude function 

Fig. 9 shows a spectrogram of the original heart sound. 

Fig. 10 shows a spectrogram of the sinusoidal coded heart sound, 
30 Fig. 1 1 shows a tracking pattern for the sinusoidal coded heart sound, 

Fig. 12 shows a Time-Frequency scaling procedure, 

Fig. 13 shows a spectrogram of a frequency shifted heart sound (p=l,g^l,r=100Hz), 
Fig. 14 shows a spectrogram of a frequency scaled heart sound (p=l,^2,r=0). 
Fig. 15 shows a spectrogram of a non-scaled sinusoidal coded heart sound 
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(p=l,^l,r=0). 

Fig. 16 shows a spectrogram of a joint time-frequency scaled heart sound (p=2,qf=2, 
r=0), 

Fig. 17 shows a spectrogram of a time scaled heart sound (p=2,q=l,r=0), and 
5 Fig. 18 shows a spectrogram of the original heart sound (p=l,q=l,i=0), but using a 
different time axis. 

The mathematical foundation for this procedure is described in the following 
paragraphs. 

10 

The sinusoidal model was originally developed for speech applications by McAiilay 
and Quantieri and a thorough description can be found in Robert J. McAulay and 
Thomas F. Quatieri, "Speech Analysis/Synthesis Based On A Sinusoidal Representa- 
tion", IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP- 

15 34, NO. 4, August 1986, page 744-754, and Thomas F. Quatieri and Robert J. 

McAulay, "Speech Transformations Based On A Sinusoidal Representation", IEEE 
Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP-34, NO. 6, 
December 1986, page 1449-1464. Because the spectral structure of a speech signal 
and a heart signal are very different (a heart beat does not have pitch like voiced 

20 speech, nor a broadband structure like unvoiced speech), the original sinusoidal 
model will not be able to model a heart beat properly. A sinusoidal model specific 
for heart sounds is developed by exploring the spectral characteristics of the heart 
beat, when tracking spectral components in the signal. 

25 The sinusoidal model is based upon sinus functions. The synthesis part of the model 
consists in adding .a group of sinusoidal fiinctions: 

Eq. (1): 

Lit) 

y(t) = ^ >4^(0cos(%(0) 

The number of functions L(t) will change as a function of time, depending on the 
30 complexity of the signal. Each sinusoidal function (frequency track) is controlled by 
the total phase function Q(0> and the amplitude A(t). The main objectives of the 
analysis part of the model is to establish the content of these two functions. 
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A fundamental relation for the sinusoidal model is the relationship between total 
phase and the instantaneous frequency: 



Eq. (2): 

0(0 = jfiiO [rad/s] 
t 

a(t) = J (a(x)dx + (p^ff^^t [rad] 

-00 

5 Substituting the expression for the total phase into eq. (1), the final expression for the 
sinusoidal model is obtained: 



Eq. (3): 

yit) = ^^^(Ocos 

Jk= 1 



The structure of a traditional sinusoidal model is illustrated in Fig. 1. The input 
10 signal x(t) is obtained from a sensor for heart sounds placed in contact with the body 
imder examination. The output signal « made available to the medically trained 
person via e.g. headphones after suitable amplification. Furthermore it may be used 
for a visual display and comparison of frequency spectra on a display screen.* 



The model may be divided into three parts: an Analysis part where the signal 
15 parameters are measured, an Interpolation and transformation part, where Si(f) and 
A(i) are interpolated from one window position to the next window position, and a 
Synthesis part where the transformed heart signaHs-reconstructedr 



The apparatus, the principle of which is shown in Fig. 1, essentially comprises means 
20 for windowing the time function, short time Fourier spectrum analysis means 1, 
means for searching and classifying spectral peaks 2, 3, means 4 for comparing 
phases of signals corresponding to said spectral peaks, means for controlling the 
phases of sine generators 6 providing signals corresponding to said spectral peaks, 
means S, 7 for controlling the amplitudes of said sine generators, and means 8 for 
2S sunmiing the signals of said controlled sine generators in order to obtain a 
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synthesized and essentially noise free output signal representative of said time 
frinction. The content of each block of Fig, 1 will be elaborated more in the following 
sections. 

5 The first step of the Analysis is to measure the frequency content of the signal. As 
the frequency content is changing as a function of time, a Short Time Fourier 
Transform (STFT) is used to track the frequency content. The STFT segments the 
signal into overlapping segments (75%), and each segment is Fourier transformed 
using a Discrete Fourier Transform (DFT). The segmentation procedure is illustrated 
10 in Fig. 2. 

When a signal segment has been Fourier transformed, the most dominant frequency 
component in the spectrum has to be found, and used as local ""check points'" for a 
frequency track over time. The original sinusoidal model finds every peak in the 

15 spectrum, and associates a sinusoidal generator to each peak. In the case of a voiced 
speech signal, the signal is periodic and the corresponding spectrum will be discrete 
(pitch). As the spectrum is discrete, the known peak search procedure will find a 
large number of peaks. However, a heart beat is non-periodic and has a nature like a 
Gaussian fianction modulated by a low-frequency sinusoid. The corresponding 

20 spectrum is also Gaussian, and the known peak search procedure will model this 
heart beat by only one sinusoidal function, where the amplitude is equal to the peak 
level. This will not model the heart beat correctly, compared to the remaining 
spectrum. However, if the amplitude of the sinusoidal function instead is set equal to 
the area below the curve defining the "hill" corresponding to the peak, the heart beat 

25 will be properly modelled. The peak searching procedure is illustrated in Fig. 3. 

First the index for the maximum peak is found, Fig. 3a. Secondly, the area of the 
"hill" is estimated. Fig. 3b, and the complete "hill" is deleted from the spectrum. Fig. 
3 c. The procedure is repeated with the modified spectrum, until a maximum number 
30 of peaks has been found, or the peak level is below a pre-defined threshold. For 
certain heart sounds as many as 50 peaks corresponding to 50 sinusoids may be 
found. 
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The result from the pecJc search is a table of peaks for each position of the window 
function. The table corresponding to window position m will be termed spec m. The 
distance between each window position is T seconds. This is illustrated in Fig. 4. 

5 The peaks in spec m now have to be matched together with peaks in spec m+I , in 
order to construct a frequency track that the total phase function £2(0 must follow. 
Many different solutions for this matching problem have been proposed in the 
literature, as no optimal solution exits. The matching strategy used in the present 
invention is based on the following three rules: 

10 - high energy peaks are matched first 

- tracks are not allowed to cross. 

- the frequency differences between two matched peaks must be below certain 
limits. 

15 First, the procedure finds the maximum peak in spec m. This peak has to be matched 
with a peak in jpec m+i. For this purpose, the procedure searches i^ec m+1 for the 
maximum peak within lower and upper limits. These limits are controlled by two 
factors. The distance to the nearest peak in spec m and the maximum allowed 
frequency change during T seconds. Before a match is accepted, it must be verified, 

20 i.e. it must not cross a previously found match. When the final match is verified, the 
two peaks are no longer visible to the matching procedure. Unmatched peaks in spec 
m are declared dead, and unmatched peaks in ^ec m+J are declared bom . The 
procedure is illustrated in Fig. 5, and an example of the matching process is 
illustrated in Fig. 6. 

25 

The synthesis will be designed for T seconds at a tune, that is from window position 
m to window position m+1, and this T second reconstructed signal will be termed 
frame m, with a frame length T. The number of bom, continued, and dead tracks in 
frame m a^II be termed £l(iw), L2(m), and L3(my For each frame wi, the bom, 
30 continued and dead tracks are synthesized separately and added to form the final 
signal: 



Eq. (4): 
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Ao=^r(o+;'2(o+;'3(o 

Xl(m) L2im) L3(im) 



Each track k in frame m has to be synthesised by a sinusoidal function 
i4jfe(t)cos(nit(/). The track is up to now only specified by the start and end points. This 
5 is illustrated in Fig. 7. 

A'^k is the amplitude of the Fourier coefficients, co^k is the frequency of the Fourier 
coefficient and is the phase offset for the current track. The phase offset is 
increased by an average phase step dependent on the start and stop frequency of the 
10 track (see Fig. 8). The original model is based on the phase of the Fourier 

coefficients. The amplitude is interpolated by ^4;^/), and the instantaneous phase is 
interpolated by Qk(t). The interpolation functions are only constrained by the start 
and stop conditions: 



Eq. (5): 








MO) = 






= ^k 










fi^'(O) = 


< 







15 

The amplitude function is not critical, so a linear interpolation function will be used 
for Au(i). However, the instantaneous phase function Clkif) is more critical because 
any discontinuity in the phase progress will deteriorate the perceptual quality of the 
final sound. In the following, the synthesis of bom, continued, and dead tracks will 
20 be designed separately. 

In the case of a continued track, a second order polynomial will be used to 
interpolate a smooth phase progress during a firame: 



25 



Eq. (6): 
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This is illustrated in Fig. 8. 



The coefficients are determined by inserting the conditions from Eq. (5) into Eq. (6): 

5 

Eq.(7): 

n(0) = ai = (pk 

n'(o) = ^2 = 

m m + 1 1 m 

The linear interpolation of the amplitude fimction is achieved by: 



Eq. (8): 

- Y — 

10 

The final equations for the synthesis of a continued track k in fi-ame m is given by: 

Continued track 
Eq. (9): 

%(0 ^ (Pjt-^^^jt — 2T 

15 

Two special cases remain to be considered. When a track is declared dead or hom^ 
we do not have both start and end conditions for them. 



In the case of a bom track, an instantaneous phase function will be initiated fi-om the 
20 begining of frame k, with^'"jk = 0, d^u = ©V \ = 0, and the second order 
instantaneous phase function is reduced to a linear function: 
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Bom track: 
Eq. (10): 

When a track is declared dead, the track is terminated at the end of frame m using 
^Wj^*^ r= 0, ^ = <x>\ and <pV ^ = 9""^ and the second order instantaneous 
S phase function is again reduced to a linear function: 



Dead track: 
Eq. (11): 



10 In order to explwn more fully the sinusoidal model according to the invention it is 
prepared for implementation in a mathematical simulation program, such as MatLab 
(TM), and the interpolation fiinctions are hence converted to discrete index. The 
discrete time version is obtfuned using the following substitutions: 



15 Eq. (12): 



t=n'A 

m 

m _ 2n 



These substitutions are now inserted into the main equations for the sinusoidal 
model. 



20 Eq. (13): 
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The time index w is running from 0 to N. The discrete amplitude function is given by: 
Eq. (14): 

[«] = • A) = ^jfc -« • A . — jyT^ ^* -« 

The final equations for a discrete-time continued track is: 

5 

Discrete continued track: 
Eq. (15): 

Ak [n] ^ A^-n — 

The discrete function for a bom and dead track can likewise be found following the 
10 same approach: 



Discrete bom track: 
Eq. (16): 

ni + l 

= /I • A • r ^'^k ' T? — 



15 DisCTete dead track: 



Eq. (17): 
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The result of coding a heart sound with severe systolic murmurs is illustrated in Figs. 
9, 10, and 1 1 . Fig, 1 1 shows the track pattern that is used to reconstruct the heart 
soimd and clearly illustrates how tracks are bom, are continued, and die, according to 
5 the spectral content of the original heart sound, and in this respect Fig. 1 1 is a 
practical demonstration of the principle shown in Fig, 6. 

Using the novel sinusoidal model described above, we are now in a position to make 
modifications to the original signal, which would normally be very difficult to 
10 obtain. As the model contains a very precise description of how the frequency 

content of the signal is changing over time, it provides a good foundation for making 
signal transformations like time stretching (without changing the frequency 
structure), frequency stretching (without changing the time progress), and joint time 
frequency stretching. 

15 

The sinusoidal model may be extended to a time-frequency scaled sinusoidal model 
with small modifications. The joint time frequency scaling is obtained by mapping 
track information at time t and frequency <d in the original model to a new time 
position r = p • ^, where p is the time scaling factor, and new frequency position 
20 (xT^q-t^ where q is the frequency scaling frictor. A pure frequency scaled model, 
may be obtained by setting = 1, and a pure time scaling model may be obtain by 
setting q = Ir The procedure is illustrated in Fig; 11. 

The time scaling is implemented in the sinusoidal model by changing the frame T to 
p' Tin the interpolation functions Q(0 and A(ty Frequency scaling can be 
25 implemented by shifting the peaks found by the peak searching procedure from co to 
q • ©: 



T ^ pT 
©- = ^ • CO 
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When q • oo crosses half the sampling frequency, the track is removed. Every time a 
bom or continued track has been synthesized, the phase ofifset (^^jt ^ has to be 
S updated to: 

Eq. (19): 

Other frequency modifications will also be possible. Heterodyne modification is 
obtained by shifting peaks from © to oo + r, where r is the amount of the spectrum 
10 movement. 

The time-scaled sinusoidal model may form the foundation for a method for auto- 
scaling a heart sound signal, which scales arbitrary heart sounds to 60 beats per 
minute (bpm). This autoscaling capability is obtained using a scale factor/?: 
15 Eq. (20): 



h 

^ ^ 60 

where h is the heart rate. Because of the sample rate in the Fourier transforms there is 
no problem in following a changing heart rate dynamically. The usefulness lies in the 
20 fact that a stressed heart will display a change in the spectral content of the heart 
sounds and focussing on tWs phenomenon may be simplified if the simultaneous 
increase in heart rate may be disregarded. 

In Figs. 13 to 18, the results of applying the above described procedure are shown 
25 when used on a heart sound. 

One particular use of the procedure relates to the transmission of very low frequency 
signals over a channel with a limited bandwidth centred much higher than than the 
low frequency signals, such as heart sound signals to be transmitted over a traditional 
30 telephone line with the bandwidth 300-3400 Hz. Traditionally this could be obtained 
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by heterodyning only, in which harmonic relationships between partials are 
destroyed. According to the present invention scaling to the frequency range of 
interest will preserve such harmonic relationships. 



5 The basic parameters that apply to the procedure described above when used with 
arbitrary heart sounds may obviously be chosen by the skilled person according to 
specific needs. However, the examples reproduced here are based on the following 
parameters: 



10 Sampling frequency: IkHz 

Window length: 64 samples 

Window shift: 4 samples 

Window type: Gaussian 

Max. number of sines: 50 

15 Max. frequency shift of tracks between neighbouring windows: 80 Hz 



It will be understood that once the signal has been converted to digital representation 
of data, its manipulation may take place in dedicated processors, RISC processors or 
general purpose computers, the outcome of the manipulation being solely dependent 

20 on the instmctions performed on the data under the control of the program written for 
the processor in order to obtain the ftinction. The physical location of the data at any 
one instant (i.e. in varying degrees of processing) may or may not be related to a 
particular block in the block diagram, but the representation of the invention in the 
form of interconnected ftinctional blocks provides the skilled person with sufficient 

25 information to obtmn the advantages of the invention. 



The foregoing description of the specific embodiments will so ftiUy reveal the 
general nature of the present invention that others skilled in the art can, by applying 
current knowledge, readily modify or adapt for various applications such specific 
30 embodiments without undue experimentation and without departing from the generic 
concept, and therefore, such adaptations and modifications should and are intended 
to be comprehended within the meaning and range of equivalents of the disclosed 
embodiments. It is to be understood that the phraseology or terminology employed 
herein is for the purpose of description and not of limitation. The means, materials. 
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and steps for carrying out various disclosed fiinctions may take a variety of forms 
without departing from the invention. 



Thus, the expressions "means to ... " and "means for or any method step 
5 language, as may be found in the specification above and/or in the claims below, 
followed by a functional statement, are intended to define and cover whatever 
structural, physical, chemical, or electrical element or structure, or whatever method 
step, which may now or in the future exist which carries out the recited functions, 
whether or not precisely equivalent to the embodiment or embodiments disclosed in 
10 the specification above, i.e., other means or steps for carrying out the same function 
can be used; and it is intended that such expressions be given their broadest 
interpretation. 
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1 . A method for arbitraiy time and/or frequency scaling of 
phonocardiographic signals obtained from a transducer and subjected to digital 
spectral analysis, characterised in that the signals are converted on a 
running basis to a sinusoidal model representation, that a time and/or frequency axis 
scaling is defined and used to control the amplitudes and phases of the sinusoids, 
which are subsequently added to create a time and/or frequency scaled representation 
of the phonocardiographic signal. 

2. A method according to claim 1, characterised in comprising the 
following steps: 

- obtaining the frequency content of the signal by applying a Short Time Fourier 
Transform with overlapping segments, a Discrete Fourier Transform being 
performed on each segment, 

- performing a frequency peak search on each segment by consecutive removal of 
the spectral components having the highest energy content, identifying each peak 
by its frequency value 

- repeating the peak search until a maximum number of peaks have been identified 
or until the energy content of the last peak is below a preset minimum, 

- establishing a segment-by-segment map of spectral peaks, said peaks forming a 
track over time 

- optionally subjecting the frequency values of the spectral peaks to a multiplication 

- defining a synthesis frame [frame m] of time, based on said segments 

- optionally subjecting each frame to a multiplication of the time scale 

. - adjusting the phase of sine generators centered on the frequencies of the.tracks 

- adjusting the amplitudes of said sine generators 

- summing the outputs of all sine generators active at any one instant for a given 
frame length T. 
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3. A method according to cl€um 1, characterised in including the 
following steps: 

- obtaining the frequency content of the signal by applying a Short Time Fourier 
Transform with overlapping segments, a Discrete Fourier Transform being 

5 performed on each segment, 

- performing a frequency peak search on each segment by repeated removal of the 
highest spectral "hills" identifying each peak by its frequency value 

- repeating the peak search until a maximum number of peaks have been identified 
or until the peak level of the last peak is below a preset minimum, 

10 - establishing a segment-by-segment map of spectral peaks, said peaks forming a 
track over time 

- adjusting the phase of sine generators centered on the frequencies represented by 
the tracks 

- summing the outputs of all sine generators active at any one instant for a given 
15 frame length T 

- creating a continuous output signal by joining consecutive frames. 

4. A method according to claim 2 for scaling a phonocardiographic signal on 
the frequency axis by a desired factor, characterised in that the spectral 

20 peaks are multiplied by a frictor q. 

5. A method according to claim 2 for scaling a phonocardiographic signal on 
the time axis by a desired factor, characterised i n that the frame length is 
multiplied by a factor/?. 
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6. A method-according-to claim 2 fonautoscaling a phonocardiographic, 
signal on the time axis, characterised i n that the scaling factor/? is set such 
that the frame length [T] multiplying factor is equal to the heart rate divided by 60. 
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7. A method according to claim 2 or 3, characterised 
number of sine generators is equal to or below 50. 



in that the 
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8. An apparatus for performing the method of any of the above claims on a 
phonocardiographic time function, 

characterised in that it comprises means for windowing the time function 
5 [xft)], short time Fourier spectrum analysis means [1], means for searching and 
classifying spectral peaks [2, 3], means [4] for comparing phases of signals 
corresponding to said spectral peaks, means for controlling the phases of sine 
generators [6] providing signals corresponding to said spectral peaks, means [5, 7] 
for controlling the amplitudes of said sine generators [6], and means [8] for summing 
10 the signals of said controlled sine generators in order to obtain a synthesized and 
essentially noise free output signal {yftj] representative of said time function. 
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